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Abstract In this study, the linear, third-order nonlin- 
ear and total absorption coefficients (ACs) of multi- 
layered quantum dot (MLQD) and multi-layered quan- 
tum anti-dot (MLQAD) with a hydrogenic impurity are 
calculated. The analytical and numerical solutions of 
Schrodinger equation for both MLQD and MLQAD, 
within the effective mass approximation and dielectric 
continuum model, are obtained. As our numerical re- 
sults indicate, an increase in the optical intensity changes 
the total AC considerably, but the intensity range that 
leads to these changes is different for MLQAD and 
MLQD. It is observed that by changing the incident 
photon energy, the AC curves corresponding to MLQAD 
and MLQD are of different shapes and behaviors. The 
peak heights of AC curves corresponding to MLQAD 
are strongly affected by changing the core antidot ra- 
dius and the shell thickness values, however in these 
cases no considerable changes are observed in peak heights 
of MLQD. Furthermore, in contrast to MLQAD, the 
photon energies corresponding to total AC peaks of 
MLQD are more affected by changing the confining po- 
tentials (CPs). 

keywords: Multi-layered quantum dots and anti-dots; 
Optical intensity; Linear and nonlinear absorbtion co- 
efficients; Confining potential. 



1 Introduction 

In the recent years, nano-structured semiconductors have 
attracted much attention, due to a wide range of their 
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applications in electronic and optoelectronic devices. 
The most known spherical nano-systems are single-layered 
quantum dots (QDs), multi-layered QDs and quantum 
anti-dots (QADs). Such systems lead to foundation of 
discrete and not degenerate energy levels like atoms [TJ 
[2ll3j|4l[5j|6l[7] . The amount of energy is dependent on the 
material, the shape, the size, the confinement poten- 
tial and the surrounding matrix of QDs. The nonlinear 
optical properties associated with optical absorption 
of mentioned nano-structures are known to be much 
stronger than bulk material, due to the quantum con- 
finement effect [gUMTOTITTllTM rS]. These particular prop- 
erties and precise engineering enable us to construct the 
optoelectronic devices, such as infrared photo detectors 
and high speed electro-optical modulators fT4lll~5lH6l 
H7] . Due to relevance of QDs to several technological 
applications, their linear and nonlinear optical proper- 
ties have been investigated both theoretically and ex- 
perimentally by many authors [TgllT*Mfll2Tll221 |%P51 



The optical properties of QDs was considered by 
Efros 27 for the first time. He studied the light ab- 
sorption coefficient in a spherical QD with infinitely 
height walls. In ref [TS], the authors have calculated not 
only the linear and nonlinear absorption coefficients but 
also the refractive index changes in a three-dimensional 
Cartesian coordinate quantum box with finite confin- 
ing potential barrier height. Ruihao Wei and Wenfang 
Xie in |19j have used the effective mass approximation 
and the perturbation theory for studying the linear and 
nonlinear optical properties of a hydrogenic impurity 
confined in a disk-shaped QD with a parabolic poten- 
tial in the presence of an electric field. Their results 
show that the optical properties are strongly affected 
by a confining strength and applied electric field inten- 
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sity. In ref [3T] , the nonlinear optical properties and the 
refractive index changes in a two-dimension (electron- 
hole system) have been investigated theoretically. The 
linear, nonlinear and total refractive index changes and 
absorption coefficients for a transitions in a spherical 
QD with parabolic potential have recently been stud- 
ied in ref [23]. Their result, expressed in several allowed 
transitions Is — lp, lp — Id and Id — If, show that the 
transition between orbitals with big I (orbital quantum 
number) values move to lower incident photon energy 
region in the presence of parabolic potential term. The 
effect of constant effective mass and position-dependent 
effective mass on the optical linear and nonlinear ab- 
sorbtion coefficient have been calculated by Khordad 
in ref [35]. 

For many semiconductor quantum heterostructures, 
such as GaAs/Gai-^Al^As, the polarization and image 
charge effects can be significant in the multi-layered sys- 
tem if there is a large dielectric discontinuity between 
the dot and the surrounding medium. However, this is 
not the case for the GaAs/Gai_ x ALjAs quantum sys- 
tem [30] , therefore these effects can be ignored safely in 
our calculation. Furthermore, for the sake of generality, 
the difference between the electronic effective mass in 
the dot (antidot) core, the shell and the bulk materi- 
als have been ignored (i.e., mi = = = m ). 
The energies are measured in meV, the effective Ryd- 
berg, Ry = me 4 /2h 2 (Aire) 2 , and distances are expressed 
in ao = An eh 2 / (me 2 ), for instance, in the particular 
case of GaAs-based semiconductors, m = 0.067me , 
and e = 13.18eo- Thus, for a GaAs host, the effective 
Rydberg is numerically Ry = 5.2meV, and the effective 
Bohr radius, ao = WAnm. 

The rest of this paper contains the following sections. 
In sec 2 we solve the Schrodinger equation for both 
MLQAD and MLQD with a hydrogenic impurity in the 
center. The solutions in each case have been expressed 
in the Whittaker and hyper geometrical functions. In 
sec 3 the linear, nonlinear and total ACs are plotted 
for various conditions as the function of incident pho- 
ton energy. Finally, our conclusions are presented in the 
last section. 



2 Theory and formulation 

The multi-layered spherical nano-systems can be made 
by binding of GaAs, Gai^Al^As and Gai-yAl^As ma- 
terials. The type of these adjacent material connec- 
tions next to each other leads to the so-called MLQD 
and MLQAD as shown in Fig. [T] More precisely, a 
MLQD (MLQAD) consists of a spherical core made of 



GaAs (Gai_ y AlyAs) surrounded by a spherical shell of 
Gai^AlzAs (Gai-zALjAs), embedded in the bulk of 
G&i- y Al y As (GaAs). In other words, the MLQAD is 
made whenever the core material and the bulk mate- 
rial in MLQD change places. We denote a for the core 
radius and b for the total dot (core plus shell) radius; 
therefore, b — a is the thickness of the shell. 

Within the framework of the effective mass approx- 
imation and in the spherical coordinate system, the 
Hamiltonian for the MLQD (MLQAD) with hydrogenic 
impurity in center is given by 

h 2 e 2 

H =-7T^V 2 --^ + V{r) (I) 

where m(r) and e are the electronic effective mass and 
dielectric constant in the semiconductor medium. m(r) 
is represented as follows 

!m\ r < a 

m-i a < r < b (2) 

m 3 r > b 

V(r) is the CP and arises due to a mismatch between 
the electronic affinities of the two regions. The barrier 
height V\ (V2) arises in layer Gai-^Al^ As (Gai-yAl^As) 
when it binds to the GaAs (Gai-^ALWs). According to 
fig- HI the potential energy V(r) in Eq. Q] for MLQD is 



(V r <a 

V(r) = I V x a<r<b (3) 

{ V 2 r>b 

and for MLQAD is 

{V 2 r <a 

Vi a<r<b (4) 

Vb r > b 



the value of Vb is always zero (i.e., the GaAs material 
have no CP) but writing Vb helps us to be able to write 
next equations in the compact form as you will see. 

Since the potential is spherically symmetric, we are 
going to use the spherical coordinates. In the three di- 
mensional spherical coordinate, the Schrodinger eigen- 
value equation is Htp(r, Q) = Eip(r, fl), where the eigen- 
functions are of the separable form ip(r, Q) — R(r)Y(f2). 
For the r-dependent (spherically symmetric) potentials, 
the angular dependence of the wave function, Y(Q), is 
given by the spherical harmonics 28 . In particular, the 
angular part is unaffected by the core and total radius, 
thus we ignore this part. The radial part is affected by 
not only the core and total radius but also the dot or 
antidot CPs, and hence, we consider it in the following 
subsections. 
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H GaAs with potential V 
I Ga Al As with potential V, 
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Fig. 1 A schematic view of the MLQD and MLQAD structures 



2.1 MLQAD 

Using the method of separation of variables in Schrodinger 
equation (1), the radial eigenvalue equation in the spher- 
ical coordinate in this case can be expressed as 



% 2 d 2 2d_ 
2m dr 2 r dr 

+ {E-Vi 



!S+V<><r) 
-)JJ«(r)=0, 



4ner'~ ^ (5) 
where i = 2, 1 and are corresponding to r < a, a < 
r < b and r > b respectively. Just as we can see in [7] , we 
can solve this equation as follows. Using the following 
definitions 



-8m(E- Vi, 



our, 



A, 



2me 



(6) 



and the change of variable as W(£i) — ^R^(r), then 
the above Schrodinger equation can be rewritten as fol- 
lowing Whittaker equation 



d 2 W&) 



6 



1)2 

11 



W{& = 0, (7) 



The general solution to this equation is the Whittaker 
functions [29] 



(8) 



where Cu and C 2 i are constants to be determined by 
continuity, asymptotic and normalization conditions. 
The asymptotic behavior indicates that the wave func- 
tion must be finite and zero in origin neighborhood and 
far away from origin, respectively. These conditions lead 
to C 22 = Cio = (since U x .j + i(£i) and M Ai i+ i(^) di- 
verge at the origin and infinity, respectively). Now one 
can write the wave functions in three regions as 



M A 2 ,Z+l(6) 
; J ' 



r <a, (9) 



6 



, a < r < b, 
(10) 



20" 



^ W+l(£o) 
Co = 



r > 6. (11) 



2.2 MLQD 

The radial eigenvalue equation (5) is also valid for MLQD 
but according to different definitions of V(r) for MLQD 
(see Eq. (3) and (4)), in this case we should use i = 0, 
1 and 2 for r < a, a < r < b and r > b regions, re- 
spectively. In the case of MLQD, the Eq. (5) has been 
solved in detail by Hsieh et al. in [B], so we do not di- 
rectly solve this equation and we shortly review their 
results for three different regions as follows 

(i=0) Inside the core of system (r < a) with Vb = 0, 
the positive energy bound states are possible for an 
electron that is confined inside the core dot; there- 
fore, solutions of the Schrodinger equation can be 
studied in three energy states [5] 
(A) E < 



where Cqa 

M is the first kind of confluent hypergeomet- 
ric function. The used parameters in the above 
solution are 



, 0A M(/ + 1-A ,2Z + 2,^)(12) 
is the normalization constant and 



6 



OA 



-8mE 



An = 



Aneh 



-2E 



(13) 
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(B) E > 

R {0) (^b) = C ob 



oo 



OB 



(14) 



Similar to the case of E < 0, Cos is the nor- 
malization constant. The recurrence relation A\ 
function and the other parameters are defined as 



4+i (AO = i, 



00 



4(A)) 

2 j MU(A>)-4Ua(A)) 



Cob 



(fc + Z)(fc-i-l) 



47T£?i V 2E 



(15) 
(16) 

k>l + 2, (17) 



m 

^- (18) 



(C) £ = 

R (0 \a oc r) = C 0C r ( - 1/2) Jai+iKc^), (19) 

Cqc is the normalization constant and J21+1 is 
the Bessel function, aoc is defined as 



«oc 



8me 2 



(20) 



(i=l) For a < r <b 

Using the convenient parameters (6), the Schrddinger 
equation in this region can be expressed in term of 
Whittaker equation as Eq. (7) which has a result 
just as the result of (10). In other words, the mid- 
dle layer in both MLQD and MLQAD behaves in a 
similar way. 

(i=2) For r > b 

In this case, similar to the above case, the Schrddinger 
equation can be expressed in terms of Whittaker 
equation as Eq. (7) but by using the asymptotic 
condition in infinity (the M function diverges at in- 
finity), the solution of equation can be written as 

# (2) (6) = C 2 e-^/ 2 e 2 U(l + 1-X 2 ,2l + 2, &) (21) 

where the parameters £2 and A2 are defined as (6). 
For reading with full detail refer to [B]. 

2.3 Linear and nonlinear absorption coefficients 

By using the density matrix formalism, the optical prop- 
erties for a spherical nano-system can be calculated. For 
this purpose the system under study can be excited by 
an external electromagnetic field of frequency u), 

E(t) = Ee iLUt + Ee~ luJt (22) 

If E is the perpendicular electromagnetic field along the 
z axis, the Hamiltonian of system in this case is Hq + 



ezE(t) where where Hq is the Hamiltonian of system 
without the electromagnetic field E(t). The linear and 
third-order nonlinear optical absorption coefficients of 
a spherical nano-system, within a two-level system (in a 
special case, from ground state to first allowed excited 
state), can be expressed as [T81I2T1I2B] 



cc>{lS) = lo 



a v nr 21 \M 21 



e R (E 21 -hLj) 2 + (nr 21 f 



(",-0 




Aa v nr 21 \M 21 \ A 



£ R "2e n r c' [{E 21 - hco) 2 + (7LT 21 ) S 



(23) 



(24) 



where, fj,, a v and n r are the permeability, carrier density 
and refractive index of the system respectively. Tiuj is the 
incident photon energy, E 2 \ — X/T 2 i, T%\ is the relaxing 
time between states 1 and 2, I is the optical intensity 
of incident wave and c is the speed of light in the free 
space. The remained quantities are defined as 



E 2 \ — E 2 — Ei, 
2 

£r — n r £ o, 

M21 = I < V^M^'i > I, 



(25) 



where M 2 % is an element of electric dipole moment ma- 
trix that in the spherical coordinate is j < ip 2 \er cos 0\il>i > 

I- 

Finally, the total absorption coefficient is 
a(ui,I) = a (1) (w) +a (3 \uj,I). (26) 

3 Results and discussion 

In this section, we calculate the linear, nonlinear and to- 
tal absorption coefficients for different shell thicknesses 
and CPs for both MLQAD and MLQD. The unchanged 
parameters used in our calculations are : Z21 = 0.2ps, 
g v = 3.0 x 10 22 , n r = 3.2, e R = 13.18. 
Figs. [21a) and (b) are related to the linear, nonlinear 
and total absorption coefficients as a function of inci- 
dent photon energy for MLQAD and MLQD respec- 
tively, when I — 20MW/m 2 , a — 2ao, b — a = 0.5a , 
Vi = 2Ry and V 2 = 5Ry. From Fig. [2J it is clear that 
the ACs of the MLQD and MLQAD have different per- 
formances with respect to incident photon energy in the 
same condition. In the case of comparing the total AC 
curves, Fig. [5] shows that the MLQD has an approxi- 
mately symmetric curve but MLQAD has a curve that 
abruptly increases then asymptotically goes to zero. In 
the case of the nonlinear AC curves, it is obvious that 
for this selected intensity value, the MLQD has no sig- 
nificant nonlinear AC but the MLQAD shows a signif- 
icant nonlinear AC. 
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Fig. 3 The linear, the nonlinear and the total ACs as a function of incident photon energy with the fixed parameters 
{a = 2ao,6 — a = 2.5ao, V\ = 2Ry and V2 = 5i?y} and six different intensities for (a) MLQAD that from up, the intensity 
value goes from I = 5MW/m 2 to 30MW/m 2 at fixed increments of 5MW/m 2 (b) MLQD that the range of intensity is 
100 < I < 600 with the fixed incremental steps of 100MW/m 2 . 



Figs. |3] (a) and (b) show the linear, nonlinear and 
total absorption coefficients as a function of the photon 
energy for MLQAD and MLQD respectively. The solid 
curves represent linear ACs where as in (22) are inde- 
pendent of intensity, the dashed (dotted) curves display 
the nonlinear (total) ACs curves. It is observed that, as 
the optical intensity increases, the total AC decreases 
for both MLQD and MLQAD. This is because the non- 
linear absorption (which is negative) enhances with an 
increase in intensity. In Fig. G2[a) that shows the be- 
havior of MLQAD, for total ACs curves, from up, the 
intensity value goes from / = bMW/m 2 to 30MW/m 2 
at fixed incremental steps of 5MW/m 2 . If this range of 
intensity is used for the MLQD, one can see that the 



MLQD has no significant changes in the nonlinear AC 
and therefore in the total AC (Fig. rjfb) is the special 
case of / = 20MW/m 2 that the MLQD has no signif- 
icant nonlinear AC) we used the range 100 < I < 600 
with a the fixed incremental steps of 100MW/m 2 as 
the Fig. 03h) shows it. In this case, at sufficiently high 
intensity, the nonlinear term causes a collapse at the 
center of the total AC peak that leads to the splitting 
of the AC peak into two peaks. 

In Fig. IHa), the total AC of both the MLQAD and 
MLQD are plotted when the shell thickness and CPs 
are fixed as b — a = la and (Vi = 2i?y, V2 — 5Ry) 
whiles a core radius has a four different values as a = 
2, 2.5, 3, 3.5ao- The graph shows that in lower incident 



(i 



A. R. Jafari, Y. Naimi 




photon energies, the total AC is mainly happened in 
MLQAD case and by increasing the amount of core 
radius, the total AC peak heights becomes larger and 
slightly shift toward larger incident photon energies. In 
the case of the MQAD, by increasing the core radius 
value, the peak heights corresponding to MLQD has no 
significant changes, but move toward the smaller inci- 
dent photon energies. Figs. 2{a) and (b) are plotted in 
the same condition, however their difference is in their 
V% values, where V2 — 5Ry in Fig. 21(a) and V2 = 00 in 
Fig-IUJb). Regarding these two figures, no considerable 
changes are observed. 

The total AC curves of the MLQAD and MLQD for 
four different shell thickness b — a = 1, 0.5, 0.25, 0.02cto, 
as the function of photon energy are plotted in Figs. [5j 
In Fig. EJa) the selected CPs and core radius are Vi = 
2Ry, V% = 5Ry and a = 2.0ao that by increasing the 
shell thickness, the peak of total AC curves correspond- 
ing to MLQAD moves to larger total AC values without 
considerable changes in incident photon energy, how- 
ever the curves related to MLQD shift to smaller energy 
regions without considerable changes in peak heights. 
Fig.EIb) shows the behaviour of the MLQAD and MLQD 
with the same conditions as Fig. [SJa) with new CPs 
(Vi = 2Ry, V2 = ooRy). By comparing the mutual 
corresponding curves for example dash-dot curves in 
Figs. [5] (a) and (b), it is observed that in the biggest 
shell thickness one can ignore the differences between 
the total AC curves but in the smaller shell thicknesses, 
the curves are more affected by changes in CPs. 

Generally, from Figs. @] and [5] it is found that un- 
der the same conditions, the total AC in MLQAD is 
greater than MLQD and also by increasing the shell 



thickness and core radius, the peak height of MLQAD 
is more sensitive and becomes significant larger, so we 
conclude that the utilization of MLQAD structures is 
more useful. 



4 Conclusion 

By studying the case of multi-layered spherical nano- 
systems with a hydrogenic impurity, the linear, third- 
order nonlinear and total optical ACs have been cal- 
culated. As our results indicates, in the case of the 
MLQAD the absorption curves abruptly increase then 
asymptotically go to zero but in the case of MLQD 
the absorption curves have symmetrical shapes. Fur- 
thermore, the range of variation of intensity that leads 
to significant changes in total AC curves is different 
for MLQAD and MLQD. It is observed that for the 
fixed shell thickness, by increasing the amount of core 
radius, the total AC peak heights become larger and 
slightly shift toward larger incident photon energies for 
MLQAD but the peak heights corresponding to MLQD 
have no significant changes but move toward the smaller 
incident photon energies. Also, in this case the change 
in CPs leads to no considerable changes in total AC 
curves. Moreover, for a fixed core radius, by increasing 
the shell thickness, the peak of total AC curves corre- 
sponding to MLQAD, moves to larger values without 
considerable changes in the incident photon energy but 
the curves related to MLQD shift to smaller energy re- 
gions without considerable changes in peak heights. Fi- 
nally, it is found that in the latter case, by decreasing 
the shell thicknesses, the total AC curves are consider- 
ably affected by CPs values. 
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Fig. 5 The total AC of MLQD and MLQAD as a function of incident photon energy with a = 2ao and four different shell 
thicknesses for (a) Vi = 2Ry and V2 = 5i?y (b) Vi = 2Ry and V2 = ooRy. 
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